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ABSTRACT 

We present the results of a careful analysis of associated absorption systems 
toward NGC 5548. Most of the well resolved narrow components in the 
associated system, defined by the Lya , C IV and N v profiles, show velocity 
separation similar (to within 10 km s~^) to the C IV doublet sphtting. We 
estimate the chance probability of occurrence of such pairs with velocity 
separation equal to C IV doublet sphtting to be 6 x 10~^. Thus it is more 
hkely that most of the narrow components are line-locked with C IV doublet 
splitting. This will mean that the radiative acceleration plays an important role 
in the kinematics of the absorbing clouds. We build grids of photoionization 
models and estimate the radiative acceleration due to all possible bound-bound 
transitions. We show that the clouds producing absorption have densities 
less than 10^ cm~^, and are in the outer regions of the broad emission line 
region (BLR). We note that the clouds which are hne- locked can not produce 
appreciable optical depths of O Vll and O Vlll, and hence can not be responsible 
for the observed ionized edges, in the soft X-ray. We discuss the implications 
of the presence of optically thin clouds in the outer regions of the BLR to the 
models of broad emission lines. 

Subject headings: galaxies: active - galaxies: individual (NGC 5548) - galaxies: 
seyfert - quasars: absorption lines 
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1. Introduction 

Associated absorption systems (i.e. absorption systems with -^abs— -2^em ), seen in 
~ 10% of the QSOs (and ~ 50% of low z AGNs), are a very good probe of the gas in 
the inner regions of the QSOs (AGNs). Studies of associated systems provide clue toward 
(i) understanding the chemical enrichment history (Petitjean et al. 1994; Hamann, 1997; 
Petitjean & Srianand, 1999), (ii) investigating the gas dynamics (inf all/outflow), and (iii) 
probing the geometry and velocity field in the broad emission line region (BLR) (Srianand 
& Shankaranarayanan, 1999). There are other ways of probing the gas in the centres 
of AGNs/QSOs. Broad emission hues and their correlated variabihty with the ionizing 
continuum probe the emission line gas (Peterson, 1993). Various observations suggest that 
in most of the low redshift, low luminosity AGNs the BLR consists of optically thick as 
well as optically thin clouds contributing to the total emission (Shields et al. 1995). It 
is realized that the ionization conditions in the optically thin clouds resembles that of 
clouds producing associated absorption lines (Shields et al. 1995). Ionization edges due to 
photoionized gas in the soft X-ray range (due to O Vll, O Vlll etc.,) are detected in large 
fraction of AGNs (Reynolds, 1997). These systems are known as "warm absorbers". The 
studies of "warm absorbers" probe a very highly ionized gas components along the line of 
sight. 

Modeling of different components discussed above requires a good knowledge of the 
shape of the ionizing radiation, location, chemical composition and density of the cloud. 
Within the allowed range of different parameters it is realized that one can get a unified 
model for the absorbing cloud which will produce the observed optical depths of ionized 
edges in soft X-ray and UV absorption lines (Mathur et al., 1995). High resolution 
observations of some of the nearby AGNs show that the associated absorption systems are 
made up of number of narrow components (Crenshaw et al. 1998). However the X-ray 
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spectra used to detect the "warm absorbers" are of very low resolution and it is impossible 
to establish any firm connection between the individual components in the UV absorption 
and the ionized edge seen in the X-rays. However, the important things one need to 
understand are (1) location of the UV absorbers, (2) kinematics of the clouds, (3) how 
much contribution such clouds will make to the emission line? and (4) how does the clouds 
respond to the continuum variation?. 

In this work we concentrate on the associated absorption systems observed along the 
line of sight to NGC 5548. We investigate the multi-component velocity structure in the 
associated system using the HST archival data. The absorption line profiles suggest that 
most of the narrow components in the associated systems, identified using the narrow 
unsaturated lines in Lya and C IV , show signatures of line-locking with velocity separation 
equal to that of C IV doublet splitting. We construct grids of photo-ionization models and 
estimate the radiative acceleration due to line absorption in various transitions. We discuss 
the nature of the absorbing clouds using the upper limit on the acceleration measured from 
the two epochs of observations, constraints on the ionization parameter, U, (which is the 
dimensionless ratio of number density of ionizing photons to that of particle density) and 
the requirements for the line-locking. 

We show that the number density of the clouds is less than 10^ cm~^. This with the 
allowed range in U suggests that the absorbing cloud will be situated in the outer regions 
of the BLR or in between BLR and NLR (Narrow emission line region). We also show that, 
under the framework of models considered in our analysis, these clouds can not produce 
appreciable optical depths of O Vll and O Vlll. Thus such clouds can not be responsible for 
the ionized edges seen in NGC 5548. We discuss the implications of radiatively accelerated, 
optically thin clouds in the outer regions of the BLR in the reverberation mapping studies. 
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2. Sub-components in the associated system 

The combined FOS (Faint Object Spectrograph) spectrum of NGC 5548, obtained for 
the AGN watch program (Korista et al. 1996), shows two distinct components in the Lya 
absorption at 1232A(i.e Zabs= 0.0134) and 1234A(i.e Zabs= 0.0151) (Mathur et al., 1995). 
The velocity separation between these components is ~ 500 km s~^, which is very close to 
the C IV doublet sphtting. These two components are clearly present in the C IV profile 
in the GHRS (Goddard High Resolution Spectrograph) spectrum (Mathur et al., 1999). 
However the low redshift system is weak and the profile is not well defined due to poor 
signal to noise. The very high spectral resolution (~ 46000) achieved by the STIS (Space 
Telescope Imaging Spectrograph) data allow us to perform a detail analysis of the individual 
velocity components. The GHRS and STIS spectra used in this study are retrieved from 
the HST archive and the pipeline calibrated data are processed using standard STSDAS 
tasks in IRAF. The original data are published already in the literature (Mathur etal ,1999; 
Crenshaw & Kraemer, 1999). 

We plot the continuum normalised STIS spectrum (dots), in the wavelength range of 
Lya absorption, in Fig |I[ As can be seen from the figure that the sub-components in the two 
systems noted above are clearly resolved. The components in the low-2; system are weak and 
well resolved, while the components in the high- 2; system are strong and blended. In order 
to resolve the components in the high- 2; system we make use of the C iv and N v profiles 
from the GHRS and STIS data. To illustrate the sub-components in this system we plot the 
C IV profile (GHRS), shifted to the Lya wavelength range, using the dashed curve. We fit 
the absorption lines using a Voigt profile fitting code (Khare et al., 1997). Minimum number 
of components needed to fit the blend is decided using the P^r degrees of freedom. The 
low- 2; system is resolved into 7 sub-components. As the continuum fit to the spectrum is 
not unique (due to unknown intrinsic emission line profile) and the cloud producing the 
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absorption need not cover the background source completely, the derived column densities 
from our fits are lower limits. However the measurement of observed wavelengths of the 
individual components are more reliable. Individual components in the low z system have 
velocity dispersion in the range 15-40 km s~^and a typical error in the determination of the 
central wavelength is ~ 10 km s^^ . The high-z system is resolved using 11 components. The 
centroid of each component is fixed to the value needed to fit the C iv and N v doublets 
from this system. The velocity dispersion of the individual components vary between 25 to 
70 km s^^and the error in the wavelength is of the order of 7-10 km 

It is interesting to note that most of the components which are resolved in the low-^ 
system have a corresponding components in the high- 2; system (as shown in Fig |T|) with a 
velocity splitting corresponding to the C IV doublet splitting within an error of 10 km s^^. 
Also there are three components in the high- 2; system which have corresponding components 
with a velocity splitting corresponding to the C iv doublet splitting among most of the 
resolved components. Thus the Lya profile seen in the STIS spectrum suggests the presence 
of line-locking with C iv velocity splitting. We estimate the probability of occurrence of 
8 pairs with velocity splitting equal to that of C iv doublet splitting (within 10 km s~^) 
among 18 clouds randomly distributed distributed in the velocity range 0-1500 km to 
be 6 X 10~^. This suggests that the line locking seen is most probably true. Thus assuming 
the existence of line-locking, in what follows we discuss various possible implications. 

The epoch of observations of NGC 5548 with GHRS and STIS are roughly 2 years 
apart. These observations can be used to get the bound on the acceleration of the individual 
components. The velocity plot of Lya , C IV and N v absorption lines observed at different 
epoch are given in Fig ^ Apart from a uniform shift of -0.05 A applied to all our GHRS 
data, we do not detect any shift in the wavelength of the centroid of lines between the data 
observed in August 24, 1996 and March 11, 1998. This gives a upper limit on the change 
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in the velocity of the components to be ~ 2 x 10~^ cm Note that this hmit is strictly 
true for the components in the high-z system only. In the case of low- 2; systems poor 
signal-to-noise in the expected region of C IV absorption in the STIS spectrum prevents us 
from making any comparison between the C IV profiles. Similarly the GHRS spectrum in 
the N V absorption region is noisy. However, the N v profile resembles very well with the 
Lya profile observed with the STIS. In the following sections we try to get the relationship 
for radiative acceleration and acceleration due to gravity of a optically thin (at the Lyman 
limit) photo-ionized cloud. 



Transfer of momentum of photons from the central source to the absorbing clouds 
through electron scattering, continuum ionization and resonance line absorption can be 
used to accelerate the clouds. If the clouds are highly ionized, as in the case of associated 
C IV absorption systems, the effect of continuum absorption is negligible. Electron 
scattering contributes to the acceleration only under high electron column densities, i.e., 
Ne > 10^"^ cm^. Thus in the range of ionization parameters and electron density one needs 
to explain the observed column densities of the associated systems the effect of radiative 
acceleration due to electron scattering and continuum absorption are negligible. Thus in 
this work we concentrate on the radiative acceleration due to line absorption only. 

It is the usual procedure to represent the ionization state of an absorbing cloud using 
a dimensionless ionization parameter, U. The distance of the cloud from the continuum 
source, r, can be written in terms of U as. 



where hh is the number density of hydrogen (H I -|- H ll) atoms, C the velocity of light. Q, 



3. 



Radiative Acceleration 




(1) 
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the number of ionizing photons emitted from the source is given by 

where L,^ is the luminosity at the frequent v. For a given spectral energy distribution Q can 
be estimated from the observed flux and a set of cosmological parameters. 



The gravitational acceleration of a cloud towards the central source is given by, 

\ Q J 



-Gm(^^^] (3) 



where M is the mass of the central engine. 

The total energy received by the cloud at any wavelength A per unit wavelength 
interval per unit time is, 

/ - {^) (^) 

where A is the area of the absorbing cloud. Total momentum available in the radiation per 
sec per unit wavelength interval at A is 

pitotal) = (5) 

This momentum is imported to the absorbing cloud through resonance line absorption 
and used to accelerate the cloud. The acceleration due to a single line transition is given by. 



LxAUuh 



(6) 



mQ 

where, r, is the optical depth of the hne m is the mass of the cloud and is the equivalent 
width to the line. If we assume the absorbing cloud to be a plane parallel slab then, 

gl = (J^I^y (7) 

\QN„mHXj ^ ' 
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where, uh, is the total hydrogen column density, mn is the mass of hydrogen atom and X is 
the mass fraction in terms of mass of hydrogen atom. In order to get the total acceleration 
we need to sum the gi' values due to all possible strong transitions. We assume X = 1.33 in 
all our models, and neglect the drag forces due to ambient medium. 



4. Line- locking 

Line-locking is a process in which the velocity separation between two absorption 
minima are equal to the velocity separation between some allowed resonance transitions. If 
the flow is dominated by radiative acceleration due to line absorption then the line-locking 
is achieved through a non-local coupling process in which photons resonant with one 
transition (at frequency ui) in one part of the flow (say at a location xi, and velocity Vi) 
Doppler shift into resonance with another transition (at frequency ^2) in another part of 
the flow (say at X2, and velocity V2) with 

V2^v^+'^^^C. (8) 

In such a case the total number of photons at 1^2 that can be absorbed by the out- flowing 
material at X2 is reduced due to absorption at Pi by the material at Xi. This reduces the 
net acceleration at X2. Thus the material at X2 will spend some time at an out-flow velocity 
corresponding to the velocity separation between ui and 1^2 . In order to keep this stable for 
long time two possible mechanisms are proposed. (1) Steady state locking: In the case of 
steady state locking the net acceleration drops to zero at X2. This requires that either the 
contribution by the line transition at z/2 dominates total radiative acceleration - and this is 
normally not the case - or that the total outward force is almost wholly counter balanced 
by an inward force (eg. gravity or pressure drag) such that their difference is smaller than 
the line radiation force due to the transition at ^2- In order to achieve such a case one needs 
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to fine tune the model parameters. In this case the hne-locking will occur at a definite 
ejection velocities in the rest frame of the ionizing source (i.e. line locking occurs on the 
absolute velocities). The double absorption troughs seen in C IV absorption in the case of 
few BALQSOs (Broad Absorption Line Quasi Stellar Objects) suggest such an accelerating 
mechanism is indeed important in some quasars (Weymann et al., 1991). (2)Non-steady 
state locking: In this case, the locking is assumed to be on velocity differences and not on 
absolute velocities relative to the central source. Locking occurs between two gas elements 
none of which is at constant velocity (vanishing acceleration) , it is only the relative velocity 
between the two that remains constant and equal to the doublet velocity. The fact that 
the line-locked pairs occur at different velocities in the case of associated narrow absorption 
line systems suggests that the non-steady state line locking is operating in these systems. 
In this case no critical balance between the outward and inward acceleration is required 
for locking to occur. What one needs is the reduction is acceleration at X2 due to lack of 
photons at V2 should be equal to the difference in the acceleration between the two fluid 
elements before locking. While neither of the velocities of the locking and locked gas, with 
respect to the central source, is constant at locking their velocity difference is. Also locking 
can occur in the in-falling as well as out- flowing clouds (Braun & Milgrom, 1989). 

5. Parameters of the Eissociated system in NGC 5548: 

In order to estimate the acceleration of the cloud one needs to get the estimate of 
various parameters related to the AGN and the absorbing cloud. In this section we try to 
get the allowed bounds on different parameters from different available observations. 
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5.1. Mass of the central engine: 

The mass of the central engine can be estimated assuming the width of the broad 
emission hnes are due to the gravitational influence of the central engine. The distance 
of the BLR from the central engine is obtained from the reverberation mapping studies. 
Krolik et al. (1991), based of the correlation between the hne of sight velocity dispersion of 
different emission lines and characteristic response time for each line, estimated the mass 
of the central engine to be in the range 10^ — 10^ Mq. Clavel et al. (1992) estimated the 
mass to be 3.7x10^ Mq, using the widths of the Lya line. Done & Krohk (1996), using the 
C IV emission line profile variabihty, estimated the central mass to be 2xlO^M0 (for two 
dimensional random motions) and SxIO^Mq (for three dimensional random motions) for 
the mass of the central engine in the case of NGC 5548. Other method of estimating the 
mass is using the accretion disk models and the continuum spectral energy distribution. Sun 
& Malkan (1989) need a minimum mass of 2 x 10^ Mq to reproduce the energy distribution 
of NGC 5548 in the framework of standard disk models. However when re-processing 
is added, the mass required reduces to 5.5x10"^ Mq. Thus we consider the mass to be 
10^ - 10^ Mq in our models for NGC 5548. 

5.2. Density in the absorbing clouds: 

ShuU & Sachs(1993) have studied the associated C IV absorption in NGC 5548 using 
the lUE monitoring data. A clear anti-correlation between the continuum luminosity and 
the C IV absorption line equivalent width is visible in their data. Using the 4 days (sampling 
interval) limit on the absorption line variability they estimated the electron density, 
ne > 5 X 10^ cm-^. Mathur et al. (1995) have studied the HST FOS spectra of NGC 5548 
observed as a part of the AGN watch program during 1993 with a sampling interval of 
~1 day. They did not find any change in the equivalent width of C IV absorption, within 
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measurement errors, in spite of a factor 1.5 variation in the UV continuum flux. However 
they have noted that the data is consistent with expected variabihty in the C IV equivalent 
width, within measurement uncertainties. The high resolution studies of NGC 5548, with 
GHRS and STIS, have revealed multi-component structure of the absorbers. Also the 
absorption lines do not cover the background source (BLR+continuum) completely and 
each component has different covering factor. The variability in the equivalent width might 
have some contribution from the varying covering factor (Srianand & Shankaranarayanan, 
1999). 

5.3. Spectral energy distribution on NGC 5548: 

It is known that the optical-UV continuum in NGC 5548 is variable and flattens when 
the object brightens. However there is little simultaneous observations of NGC 5548 in 
various wave bands and the change in the spectral energy distribution as a function of 
source luminosity is unknown. Dumont et al.(1998) have provided detailed discussion on 
spectral energy distribution of NGC 5548. In this study we use the continuum used by 
Dumont et al (1998). Assuming = 0.5, = 75 kms~^ Mpc~^ and redshift of NGC 5548 
to be 0.0175 we estimate the value of Q = 3 x 10^*^5"^ from the flux in the range 1 to 20 
Ryd. Note that the spectrum we reconstruct from the observations need not be the one seen 
by the absorbing clouds. As noted by Dumont et al. (1998), we should keep in mind that 
whole continuum distribution in the observed range is determined within about a factor two 
of uncertainty ( and a large factor in the Extreme UV range), owing to temporal variation. 



5.4. Ion column densities: 



Ultraviolet spectra of NGC 5548 observed with the Hopkins Ultraviolet Telescope 
(Kriss et al. 1997) show no Lyman limit at the systemic redshift. This suggest that 
N{H l) < 10^^ cm~^ in individual components. The Voigt profile fits gives lower limit on 
the N{H l) as we do not know the exact value of the covering factor. The N{H l) in the 
individual components obtained from the profile fits vary between 10^^ to 2 x 10^^ cm~^. 
From the absence of C II absorption we estimate the 2a upper limit on N{C II ) to be 
10^^cm~^. The lower limit on the C IV column density derived from different components 
in the high-z system varies between 2x10^^ to 10^^ cm . Similarly, the lower limit on the 
column density of N v is in the range 2x10-*^^ to 4x10-*^^ cm~^. It is interesting to note that 
N V column densities in all the components are higher than C IV . If the heavy element 
abundance ratios are like solar this will mean high degree of ionization. However excess N 
abundances are known in high-2; QSOs (Hamann, 1999; Petitjean & Srianand, 1999) and 
the observed ratio could be due to the excess [N/C]. Other possibility for getting the higher 
values of N(N v) could be the difference in the covering factor for N v and C IV (Petitjean 
& Srianand, 1999). ROSAT PSPC observations by Nandra et al.(1993) have confirmed 
the soft X-ray absorption by ionized gas in NGC 5548. Reynolds (1997) has obtained 
To VII — 0.24±0.04 and tq vm = 0.16±0.03 by fitting the observed ASCA spectrum using 
r = 1.88 ± 0.01 and Nh < 2.2 10^^ cm'^ The model fits to the ASCA data investigated 
by George et al.(1998) suggest that the warm absorbing material is completely covering 
the X-ray continuum emitting region. If we assume complete coverage then the estimated 
column densities are N{0 Vll) = 9.3 x 10^^ and N{0 Vlll) = 1.2 x 10^* cm"^. Done et al. 
(1995) using the ROSAT spectra obtained a lower limit on the electron density in the warm 
absorbers to be 5x10^ cm~^. Unhke UV spectrum the X-ray spectra is of lower resolution 
and it is very difficult to get the optical depth contributed by individual components, we 
detect in the UV spectra. 
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6. Results of the model calculations: 

We use the photo- ionization code cloudy (Ferland,1996) to get the ionization conditions 
in a typical absorbing cloud. We consider the absorbing cloud to be plane parallel slab 
of uniform density with chemical composition similar to that of solar composition. The 
ionizing radiation is assumed to be similar to the one used by Dumont et al(1998) and 
assume to shine on the cloud from one side. We run grids of photo-ionization models 
covering -2.0< log U <1.00 and stopped the code when the column density of neutral 
hydrogen reaches a given value. In this work we discuss the results obtained when N{H l) 
is 10^^ (model A) and 10^^ cm~^ (model B) respectively. The model results are plotted in 
Fig 0. The panels (a) and (b) give the column densities of various ionization states that are 
detected in the associated system as a function of ionization parameter. In the range of 
ionization parameter considered here, the C ll column density is well within the upper limit 
we get from the non-detection. Also the resulting N v column densities are higher than 
that of C IV . The horizontal lines in both the panels show the lower limit on the column 
densities of O vill and O vil so that the optical depth is equal to one tenth of the measured 
optical depth. Fig |^ suggests, in order for UV absorber to contribute at least 10% of the 
optical depth of Oxygen edges, that log U > —0.2 (model A) and log U > —0.6 (model B). 
If only one of the clouds produce all the observed optical depth of ionized Oxygen edges 
then log U > 0.6 (model A) and log U > —0.2 (model B). In the case of model A, such a 
cloud will not produce detectable amount of C IV and N v column densities as in the case 
of the sub-components in the low-z system (if the abundance ratios are solar). However in 
the case of model B, in a narrow range of ionization parameter, one can get a consistent 
solution for the UV as well as "warm absorbers". 

The column densities of various ionization states of H, He, C, N and O, resulting 
from our models, are stored in an array and are used to estimate the equivalent widths of 
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different transitions. We also store the total column density of hydrogen resulting from 
the models. We consider the Doppler broadening and natural broadening of lines (Voigt 
profiles) while calculating the equivalent widths. We estimate the radiative acceleration due 
to line absorption using the column densities of H I, He I, He ll and all possible ionization 
states of O, N and C. We use i/o = 75 kms-^ Mpc-\ Q = 3x10^^ s^\qo = 0.5, redshift 
of the AGN, z = 0.0175. The results presented in Fig |^ assume the velocity dispersion to 
be 25 km s~^. The acceleration due to gravity and the radiation pressure are calculated 
assuming the mass of the central engine to be 10^ Mq and lO^M©. The net acceleration 
as a function of ionization parameter are plotted for a range of densities for the two cases 
(see panels c-f Fig.^). The allowed range in the acceleration, from the observations, are 
given by the horizontal dotted lines. Positive values of the acceleration means the radiative 
acceleration dominates and the negative values represent the gravitational deceleration. 

When mass of the central engine is lO^M©, radiative acceleration dominates when the 
ionization parameter log U < —0.75 (model A) and log U < —1.45 (model B). Note that 
such a cloud will have the optical depths of O vil and O Vlll less than 0.01 and could 
not make appreciable contributions to the observed optical depths. Decreasing the mass 
of the central source will increase the upper limit on log U . Note that when we consider 
M = IO^Mq a small range is allowed in log U where the cloud can contribute one tenth of 
the observed optical depth of the Oxygen edges. However if all the Oxygen edge optical 
depths are due to single cloud then the cloud's motion will be dominated by the radiative 
acceleration due to line absorption only when M < lO^M©. Such a low values of the mass 
are not allowed by the observations. 

When we consider uh = 10^ cm~^, observational constraints on the change in velocity 
allows an extremely narrow range in U when N{H l) < 10^^ cm~'^. Also this demands a near 
critical balance between the radiative acceleration and the gravitational deceleration, which 
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requires extreme fine tuning. Also such a balance can not be maintained when the ionizing 
continuum changes as in the case of NGC 5548. Thus it is more likely that the density of 
the absorbing cloud has to be less than 10^ cm~^ when N{H l) < 10^^ cm~^. This will 
mean the distance of the cloud from the central source has to be higher than 20 It days (for 
log U < —0.2). Transfer function derived from the HST monitoring data suggest that the 
region of C IV emitting clouds streaches from 1 It-day or closer, out to 10 It-day are further 
(Done & Krolik, 1996). Thus the locations of the clouds with N(H l)< 10^^ cm~^ will be 
the outer regions of the BLR. Note that the Voigt profile fits to the components suggest 
that the N{H l) in individual components are less than 10^^ cm'"^ . However, as discussed 
before, due to unknown degree of saturation and covering factor the column density derived 
will be a lower limit only. The Lya line of the low- 2; system is well resolved and, as can be 
seen from Fig |l], the lines do not show any signature of saturation (flat bottom) and it is 
more likely that the N{H l) < 10^^ cm^^ in the individual components. Thus the location 
of such a cloud will be between the BLR and NLR. However the change in the velocity for 
these clouds are not measured due to poor S/N in the GHRS & STIS data. 

When we consider N{H l) = 10^^ cmr"^, critical balance is needed only for 
nu > 10^° cmT^. This means the distance of such a cloud from the central engine has to be 
higher than 2 It days (for log U < 0.2). Note that some of the components in the high- 2; 
system could have such large column densities. However the C iv and N v profiles of these 
components suggest that these clouds cover more than 80% of the BLR and have to be 
outside the BLR. Thus we believe A^^ < 10^ cm~^ even in these clouds. 

In all our calculations we consider the effect of the continuum radiation alone. However 
if the absorbing clouds are outside the BLR then the line photons from the emission line 
clouds could also contribute to the radiative accleration. Specially C iv and Lya emission 
lines are very strong and the excess contribution to the accleration due to these lines 
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could be appreciable. Also this effect will be larger for the high- 2; system. In order to 
keep our analysis simple we do not include the effect of the emission lines in the radiative 
acceleration. 

7. Implications of line-locking in NGC 5548 

In order for the line locking to occur the radiative acceleration should contribute 
appreciably to the net acceleration of the cloud. As discussed before the locking among 
different sub-components in absorption occures at the velocity difference corresponding to 
the velocity splitting of C iv doublet. This means the fractional acceleration contributed by 
the C IVA1550 absorption line is important. Since the locking occures at different velocities 
with respect to AGN it is clear that the case we have to consider is non-steady state line 
locking. Also if the continuum source varies it will be easy to preserve the locking under 
the non-steady state conditions. In Fig. |^ we plot the fractional acceleration contributed by 
the C IVA1550 absorption in different case under consideration. As most of the important 
transitions are saturated effect of velocity dispersion is important. In Fig ^, we plot the 
results for two different values of the velocity dispersion parameter in order to get the feel 
for the dependence of radiative acceleration on our choice of velocity dispersion. Note that 
these curves are independent of the choice of uh- In the range of ionization parameter 
where the radiative acceleration dominates C lvA1550 absorption contributes up to 4% of 
the net acceleration. The maximum contribution of the C lvA1550 line to the acceleration 
occurs when there is a near balance between the gravitational acceleration and the radiative 
acceleration. Locking at such ionization parameters are very unstable against a small 
changes in the intensity of the ionizing continuum. 

Let Oi and a2 are the net acceleration of two clouds along our line of sight (with 
0-2 > o-i)- Let fc IV be the fractional acceleration contributed by the C lvA1550 absorption 
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in cloud 2. At the locking, 



Ol = (1 - fc IV ) 0.2 



(9) 



If the ionizing continuum is same for both the clouds (i.e. cloud 1 is optically thin to 
the continuum radiation) then the difference in acceleration between two clouds can be 
achieved by changing U, uh, N{H l) and Nh- Note that the total acceleration is directly 
proportional to the produce Uuh- In order to achieve the locking the cloud which is getting 
locked (cloud 2) should be further away from the source compared to the cloud which 
produces locking (cloud 1). This demands, 



where superscripts a and b are used to represent the first and second cloud respectively. 
This condition suggests that if the ionization parameter of the locking cloud is less than or 
equal to that of the locked cloud (i.e C/" < U^) then the density of the locking cloud should 
be higher than that of the locked cloud (n^ > n^j). In order to understand the range of 
parameters one needs to get the locking we consider the following cases. 

• case i: If we assume N(H l) for locking cloud is less than or equal to that of the locked 
cloud, and if is less than then ai will be always greater than 02 in the range 
of ionization parameters relevant to this study. Thus in order for the locking to occur 



• case ii: If we assume N(H l) for locking cloud is greater than or equal to that of 
the locked cloud, then 02 will be much higher than ai. However in order to produce 
locking at the C IV velocity splitting, the difference in the acceleration between the 
clouds should be equal to the acceleration of cloud 'b' due to C lvA1550 absorption. 
This condition suggest, in the range of U relevant to this study, will be equal to 



Hi 



(10) 
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or slightly less than in order for the line-locking to occur at the C iv doublet 
splitting. 

In the case of NGC 5548 the locked clouds (low-2;) most probably have a lower H I 
column density than the locking cloud (high- 2;). As discussed above the uh of the low- 2; 
system should be equal to or less than that of the high- 2; system in order to produce 
line-locking. Though we do not have stringent limit on the change in the velocity of 
the low- 2; system, due to poor S/N, in the wavelength range of the C iv absorption, the 
above condition suggests that the density of the low- 2; system is less than 10^ cm~^. From 
Fig ^ it is clear that in the range of ionization parameter, where the radiative acceleration 
dominates, the optical depth of O vil and O Vlll in the clouds are much lower than the 
observed value. Though low values of central mass can produce solutions with radiatively 
accelerated clouds with appreciable optical depths of O vil and O vill, the fractional 
acceleration contributed by C lvA1550 in such a clouds are negligible. Thus the clouds 
which are showing line-locking with C iv doublet splitting are most probably not responsible 
for the observed ionized Oxygen edges. 



8. Discussion: 

HST observations of nearby AGNs suggest that ~ 50% of the seyfert 1 galaxies do 
show associated absorption systems (Crenshaw et al., 1998). Invariably, in all cases ^abs is 
less than , suggesting a component of out-flowing, highly ionized, optically thin gas, 
close to the centers of seyfert galaxies. The strong absorption components have cores that 
are much deeper than the continuum flux levels, indicating that the regions responsible for 
these components lie completely outside the BLR. This result is supported by the derived 
covering factors for few systems (NGC 3516, NGC 3783, NGC 4151, NGC 5548, Mrk 509) 
using the high resolution data (Crenshaw et al., 1998), and the constancy of the wavelength 
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of the absorption lines over a period of few years seen in few systems (NGC 5548 (this work) 
and NGC 4151 (Weymann et al. 1997)). Presence of signatures of hne- locking, in NGC 
5548, suggest that the radiative acceleration plays (or played) a vital role in accelerating 
the clouds. Note that even in the case of NGC 4151 some of the narrow components seem 
to have velocity sphtting close to the C IV doublet sphtting (see Fig 5. of Weymann et al. 
1997). 

The 1993 monitoring campaign on NGC 5548 (Korista et al. 1996) revealed a red/blue 
symmetry between the blue and red line core light curves of the C IV emission. However no 
such relationship is seen between the blue and red wing light curve suggesting a possible 
radial motions for emitting clouds with large velocities. Done & Krohk (1996), noted that 
the response function for the red and blue cores are nearly identical, and the red wing 
response has a strong peak at small lags but the blue wings response raises only very 
slightly toward zero time delay. This indicates, if we assume all the emitting cloud to have 
a positive response, that there is more material near the line of sight and/or at small radii 
that is traveling away from us than is traveling towards us. 

Detail studies of UV emission lines in NGC 5548 provide some indirect evidence for the 
presence of optically thin line emitting clouds. Sparke(1993) finds that gas with negative 
response is required, at small time delays, to explain the shape of the cross-correlation 
function for several lines in the case of NGC 5548. Koen(1993) finds evidence for negative 
instantaneous response of Si IV and C ll lines to the UV continuum variations.Various 
implications of the presence of such optically thin components are discussed by Shields 
et al. (1995). Observations discussed here clearly suggest that there is a component of 
optically thin gas, which are most probably present just outside the BLR, with net outward 
velocity. As the densities are high the cloud will respond to the continuum variations 
instantaneously. The C IV line emission from these clouds will have negative response to 
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the continuum variations. If the volume fiUing factor of such clouds are high then they are 
expected to contribute appreciable to the total C IV emission. The slow response towards 
the zero time lag in the case of blue wing could be due to the instantaneous negative 
response of small fraction of optically thin clouds close to our hne of sight. If these clouds 
are far away from the central continuum source and have radial outflows, like the individual 
components in the absorption systems, then they will affect only the blue wing response 
function at short time delays. However these clouds can contribute appreciable amount to 
the response function of the red wing at higher time delay. 

It is interesting to note that the response functions derived for total C IV as well as red 
wing of the C IV emission show minima (even negative response) at a delay time of ~ 20 
days (Done & Krolik, 1995). Similar result is also obtained by Wanders et al. (1995) using 
different numerical method. They interpreted this as an evidence for bi-conical structure 
in the emission line region. However this result is also consistent with our expectation in 
the case of presence of radially out-flowing optically thin gas in the outer regions of the 
BLR. Though there arc models of BLR with optically thin as well as optically thick clouds 
are discussed in the literature (Goad et al. 1995), optically thin clouds are considered very 
close to the central source. It will be an interesting and important exercise to build models 
of BLR with optically thin clouds in the outer regions as suggested by the absorption line 
studies. 



9. Summciry 

We study the multiple velocity component structure in the associated absorption 
system in NGC 5548 using HST high resolution spectra. Most of the narrow components 
in the associated systems show the signature of line-locking with velocity separation equal 
to C IV doublet splitting. Using the spectra observed at two different epochs separated by 
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~ 2 yrs we estimate the upper limit on the acceleration of individual clouds. We construct 
grids of photo-ionization models and estimate radiative acceleration due to line absorption 
by various transitions. Using the constraints on the ionization parameter, acceleration, 
covering factor and the conditions to be satisfied for the line-locking, we show that (a) the 
density of the cloud is less than 10^ cm~^, (b) the clouds are situated out side the BLR, 
(3) the clouds involved in the line-locking can not be responsible for the observed ionized 
Oxygen edges in the soft X-ray. We discuss the implications of the presence of optically thin 
clouds in the outer regions of the BLR. Lack of information on the ionization conditions of 
the clouds prevents us from making a more reahstic models of the line-locked clouds. We 
believe, future UV observations of NGC 5548 using FUSE will provide better constraints on 
the ionization parameters of the individual components, which can be used to get a more 
reahstic models. 
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Fig. 1. — Lya absorption line of the associated system toward NGC 5548. The dots and 
continuous curves represent the observed spectrum and Voigt profile fits respectively. The 
dotted curve gives the profile of the C lvA1548 line (suitably shifted to the Lya wavelength). 
Some of the components which are separated by the C iv doublet sphtting are also marked 
in the figure. 
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Fig. 2. — Profiles of Lya , C iv and N v absorption lines observed at different epoch. The 
outflow velocity is estimated assuming the redshift of NGC 5548 to be 0.0175. 
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Fig. 3. — Results of estimation of radiative acceleration using photoionization models. 
Panels (a) & (b) give the column densities of various ions as a function of ionization parameter 
for N(H l) = 10^^ and 10^^ cmr'^ respectively. The net acceleration estimated for different 
models are given in panels c(M = lO^M©, N(H l) = 10^^ cm-^), d(M = lO^M©, N(H l) = 
10^^ cm-2), e(M = lO^M©, N(H l) = lO^^ cmr'^) and f(M = lO^M©, N(H l) = 10^^ cmr'^). 
The horizontal dashed lines gives the observed hmit on the acceleration. The results are 
presented for four different value of densities. The values of log uh are marked close to the 
corresponding curves. 
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Fig. 4. — Fractional acceleration contributed by the C lvA1550 absorption as a 
function of ionization parameter. The solid and dotted curves are the results obtained 
assuming the values of velocity dispersion in the absorbing clouds to be 25 km s~^and 10 
km s~-'^respectively. 



